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We describe here a minimal theory of tight binding electrons moving on the square planar Cu lat- 
tice of the hole-doped cuprates and mixed quantum mechanically with pairs of them (Cooper pairs). 
Superconductivity occurring at the transition temperature T c is the long-range, d-wave symmetry 
phase coherence of these Cooper pairs. Fluctuations necessarily associated with incipient long-range 
superconducting order have a generic large-distance behaviour near T c . We calculate the spectral 
density of electrons coupled to such Cooper pair fluctuations and show that features observed in 
Angle Resolved Photo Emission Spectroscopy (ARPES) experiments on different cuprates above 
T c as a function of doping and temperature emerge naturally in this description. These include 
'Fermi arcs' with temperature-dependent length and an antinodal pseudogap which fills up linearly 
as the temperature increases towards the pseudogap temperature. Our results agree quantitatively 
with experiment. Below T c , the effects of nonzero superfluid density and thermal fluctuations are 
calculated and compared successfully with some recent ARPES experiments, especially the observed 
bending or deviation of the superconducting gap from the canonical d-wave form. 



I. INTRODUCTION 

High-temperature superconductivity in hole-doped 
cuprates, accompanied by a 'pseudogap phase' as well as 
other strange phenomena, continues to be an outstand- 
ing problem in condensed matter physics for a quarter 
of a century now [1-3]. Over the years ARPES [4, 5] 
has uncovered a number of unusual spectral properties 
of electrons near the Fermi energy with definite in-plane 
momenta. This low-energy electronic excitation spec- 
trum is of paramount importance for explaining the rich 
and poorly understood collection of experimental find- 
ings [1—3] from thermodynamic, transport and spectro- 
scopic measurements on cuprate superconductors. We 
show here that the spectral function of electrons with 
momentum ranging over the putative Fermi surface (re- 
covered at high temperatures above the pseudogap tem- 
perature scale [6-8] ) is strongly affected by their coupling 
to Cooper pairs. On approaching T c i.e. the tempera- 
ture at which the Cooper pair phase stiffness becomes 
nonzero, the inevitable coupling of electrons with long 
wavelength (d-wave symmetry) phase fluctuations leads 
to the observed characteristic low-energy behavior ob- 
served in ARPES experiments. The idea that Cooper 
pair phase fluctuations are important in cuprate super- 
conductivity has a long history; we mention only a few 
examples. The experimental realization that Cooper pair 
phase fluctuations are significant for a large regime of 
hole doping (below optimum doping) and temperature 
owes largely to the observation of large Nernst effect [9] 
and enhanced fluctuation diamagnetism [10] by Ong and 
coworkers. The same physics is implied in the early the- 
oretical work [11] of Emery and Kivclson. 

The cuprates [2] (e.g. La2- x Sr x Cu04) exhibit super- 
conductivity at unusually high temperatures on doping 
the parent compound, a Mott insulator, with holes (x 
per Cu site in the above case; Fig. la shows the square 



Cu lattice with lattice spacing a) . There is long-standing 
evidence from the (cosfc^a — cosk y a) dependence of the 
superconducting gap Ak on the in-plane momentum k 
of the low energy electronic excitations (observed, for 
example, in ARPES experiments [4, 5]) that the pair- 
ing involves nearest neighbours on a tight binding lat- 
tice. We therefore assume that the basic Cooper pair 
in hole-doped cuprates is the nearest-neighbour spin sin- 
glet. Superconductivity is the long-range phase coher- 
ence of these Cooper pairs. The low energy degrees of 
freedom then are (the electrons and) the complex singlet 
pair amplitudes 

= = \( a ii a rt ~ a^a 3i ) EE Vm- (1) 

Here a,i a (a\ a ) destroys (creates) an electron with spin a 
at the lattice site i, sites i and j are nearest neighbours 
and the bond between them is uniquely labelled by the 
bond-centre site m (see Fig. la). An interaction of the 
form CRe (^> m ?/>*) between nearest-neighbour bonds cen- 
tered at m and n, with positive C, leads to the observed 
d-wave symmetry superconductivity [12] below T c . The 
crossover temperature T* (x) below which the equilibrium 
value (IV'ml) of the local pair amplitude becomes signifi- 
cant is taken to be the 'pseudogap' temperature. There 
is considerable experimental evidence for this view [6- 
8], though there is also the alternative view that T*(x) 
is associated with a new long-range order, e.g. (i-density 
wave [13], time reversal symmetry breaking circulating 
currents [14], electron nematic order [15], stripes[16] etc. 

We work out here the effects of the quantum mechani- 
cal coupling between (fermionic) electrons and (bosonic) 
Cooper pairs of the same electrons. On approaching T c 
from above, collective d-wave symmetry superconduct- 
ing correlations develop among the pairs with a charac- 
teristic superconducting coherence length scale £ which 
diverges at the second-order transition temperature T c . 
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FIG. 1: Square lattice where cuprate superconduc- 
tivity primarily resides and the pair fluctuation ex- 
change process for the self energy of electrons hop- 
ping on it. a, The blue circles at {Ri} represent the Cu 
atoms and the magenta ones at {R m } (where R m = R^ = 
Ri M with /x = x,y) represent bond centres (also the location 
of oxygen atoms). The arrows denote planar spins equiva- 
lent to tpij = i'i^ = 4>m = A m exp(i# m ) with A m being the 
length of the 'spin' and 9 m the angle it makes with the x- 
axis; d-wave superconductivity translates in spin language to 
Ne'el order as shown, b, Feynman diagram[26, 27] of the pro- 
cess in which an electron (k,ioj n ) virtually becomes another 
electron (— k + 2q, — iu n + iz m ) and absorbs a Cooper pair 
(2q, iz m ) in the intermediate state. The curly line denotes 
the Cooper pair propagator D M „(2q, iz m ) (here z m — 2mn / 'ft 
is the bosonic Matsubara frequency, m being an integer). 



These correlations have a generic form at large distances 
(>> the lattice spacing a). As we show here, the effect of 
these correlations on the electrons leads, for example, to 
temperature-dependent Fermi arcs [17] and to the filling 
of the antinodal pscudogap in the manner observed [18]. 
Further, the observed long-range order (LRO) below T c 
leads to a sharp antinodal spectral feature [19, 20] re- 
lated to the nonzero superfluid density, and thermal pair 
fluctuations cause a deviation ('bending') of the inferred 
'gap' as a function of k from the expected d-wave form 
(cos k x a — cos k y a). The bending, being of thermal ori- 
gin, decreases with decreasing temperature, in agreement 
with some ARPES measurements [21, 22]. 

We describe the model used for our calculation in Sec- 
tion II, followed by a section (Section III) on the method 
of calculation and phenomenological inputs that are ex- 
tracted from various experiments. We have compared 
our results with a large number ARPES experiments, as 
mentioned in the preceding paragraph. These are re- 
ported in Section IV. There have been many previous 
studies of the effects of phase fluctuations on electrons; 
we mention a few of them in the discussion section (Sec- 
tion V) along with some concluding remarks. Our at- 
tempt here is based on a particular view of the supercon- 
ducting transition, does not assume preformed d-wave 
pairs, and a detailed comparison with a large number of 
ARPES experiments is made. The present paper, in our 
view, represents a significant development, an extension 



of the approach of Ref. [12] to explicitly include electronic 
(fermion) degrees of freedom in addition to the Cooper 
pair (bosonic) ones. It also shows that a simple mini- 
mal theory can account for for a variety of ARPES data, 
which have hitherto been analyzed separately, using phe- 
nomenological [23] or problem specific [24, 25] models. 
Appendices describe some technical details of our calcu- 
lations. 



II. MODEL 

A complete model for electrons in the square planar Cu 
lattice with relatively weak interlayer coupling and with 
interactions which strongly affect their motion is needed 
for describing low energy cuprate phenomena comprehen- 
sively. There is no such unanimously accepted descrip- 
tion yet. We use a partial, phenomenological model as 
follows. The electrons are described by the Hamiltonian, 

H = - X) *« a l<r°jV - ■''i' , >j h :r ( 2 ) 

i,j,cr <ij> 

which consists of a quantum mechanical intcrsite hopping 
term, where ty is the effective amplitude for an electron 
to hop from site i to site j, and a nearest neighbor pair at- 
traction term with strength Jy . The parameters (i.e. ty 
and Jy) of both are strongly affected by correlations. 

In the case of strong local repulsion such as large Mott- 
Hubbard U (» t) [28], one uses the same form as in its 
absence, but with renormalized Hamiltonian parameters 
which are assumed to describe the entire effect of U (at 
least at low energies). This implies (for example with 
renormalized hopping fy replacing ty) that there are 
good mobile quasiparticles, albeit with renormalized dis- 
persion. The well known single-site Gutzwillcr rcnormal- 
ization factor [28-30], for example, projects out doubly 
occupied sites completely and is a good approximation 
for U — > oo. It assumes, for small hole density x, a multi- 
plicative factor g t = 2x/ (1+x) for ty and g s = 4/(1 + x) 2 
for Jy (i.e. fy = gttij and Jy = g s Jij)- The bare hop- 
ping amplitude ty involves the nearest neighbor (t), next- 
nearest neighbor (t') and further neighbor (f") hopping 
terms. In our calculations, we use the above mentioned 
homogeneous Gutzwillcr approximation with standard 
values for ty (f = 300 meV, t'/t = -1/4 and t" = 0, see 
e.g. [31]). The value of Jy does not appear explicitly in 
the calculation, as we elucidate later. The Jij term is the 
well known superexchange (e.g. Jy oc tfj/U in single site 
Hubbard model [28]) for large U. We have written the 
conventional JySj.Sj term as a pair attraction — Jy&y&y 
[32, 33] using the identity (Si.Sj— fiifij/ty = — 6y6y valid 
for spin-i particles (S, and hi are the spin and number 
operators at site i, respectively). There is thus strong at- 
tractive nearest-neighbor spin singlet pairing term in the 
Hamiltonian, given that the antifcrromagnetic coupling 
Jij is known to be large (~ 1500 K) [34] for undoped 
cupratcs (x = 0). 
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The above identity between nearest-neighbor AF 
Heisenberg spin interaction and spin singlet Cooper pair 
attraction means that exact solution of Eq. (2) with cither 
term would lead to the same result. Because of the fact 
that antiferromagnetic LRO disappears for surprisingly 
small hole doping x and is replaced by superconducting 
order, and that the holes are quite mobile, we use the 
singlet bond pair form in Eq.(2) as being natural and ac- 
curate for good approximations, e.g. mean field theory, 
and assume a homogeneous system as will naturally arise 
for mobile holes. 

Very generally, the bond-pair self interaction term 
in Eq.(2) can be written via the exact Hubbard- 
Stratonovich transformation [35] as a time (and space) 
dependent bond pair potential acting on electrons and 
characterized by a field tp m {T) with a Gaussian probabil- 
ity distribution [here r is the imaginary time; < r < 
j3 = l/(fcsT)]. The saddle point of the resulting action 
in the static limit gives rise to the conventional mean 
field approximation in which the second term in Eq.(2) 
is written as 

- ^Mi)K +b\ 3 {b l0 ) - {bl){b l3 )) (3) 

and the average ip m = tpij oc J%j(bij) is determined self- 
consistently (mean field theory). 

The effective Hamiltonian we use (see Section III) is 
of the form of Eq.(2) with the second term in it replaced 
by Eq.(3). This describes two coupled fluids, namely 
a fermionic fluid and a bosonic fluid, represented re- 
spectively by the on-site electron field a\ a and the bond 
Cooper pair field ipy = ip m . The properties of ip m needed 
in our calculation are its mean value (ipm) (nonzero below 
T c ) and the fluctuation part of the correlation function 
(ipmijj*) (whose universal form for large |R m — R„| near 
T c is what we use). These arise from inter-site interac- 
tions of the Cooper pairs tp m ; our results are independent 
of the details of these interactions. A nearest neighbor 
interaction of the form C J2 <mn > (V'mV'n + ipm^n)i with 
C > 0, has been mentioned earlier as a natural possibility 
[12]. It may arise [12] with C oc x in a strong correlation 
picture due to diagonal hole hopping t' . 

For a translationally invariant system described by 
Eq.(2), the electron Green's function satisfies the Dyson 
equation [26, 27] 

G-\k,iuj n ) = (G°)- 1 (k,^„)-S(k,i w „), (4) 

where E(k, iuj n ) (oj n = (2n + l)ir/(3 is the fermionic 
Matsubara frequency with n as an integer) is the ir- 
reducible self energy, originating from the coupling be- 
tween bond pairs and electrons with bare propagator 
G°(k,MJ n ) = (iw n - Ck) _1 where £ k = e k - (j, and e k 
is the Fourier transform of the hopping t%j = (/z is 
the chemical potential). 

We use a well-known bosonic fluctuation exchange ap- 
proximation for E(k, iu) n ) which captures the leading ef- 
fect of this coupling beyond the 'Hartree' approximation 
and is shown diagrammatically in Fig. lb. This diagram 



describes the exchange of a Cooper pair fluctuation by 
an electron. In common with general practice, we find E 
(and thence G) by inserting G° instead of G in the expres- 
sion for it [see Eq.(5) below]. This is known to be gen- 
erally quite accurate [26], e.g. for the coupled electron- 
phonon system. 

Close to T c , the temporal decay of long- wavelength 
fluctuations is specially slow (dynamical critical slowing 
down [36]) so that they can be regarded as quasistatic 
[12] (i.e. the characteristic frequency scale uj « k B T). 
Though the decay is 'slow', we assume, as is generally 
done, that the system is homogeneous. There are an- 
nealing processes which make it so on the experimental 
time scale. (There may be intrinsic as well as extrinsic 
sources of static quenched disorder, e.g., due to the very 
process of doping itself [37]). Even their effect can be 
included by approximate configuration averaging. Our 
approach is similar to that for static critical phenomena 
where also a homogeneous system is used and all the 
fluctuations are thermal. The self energy E can then be 
expressed as 

E(k,iw„) = 

E G°(-k + 2q,^„)^(2q)/ Al (k,q)/,(k,q) 

(5) 

where TV is the total number of Cu sites on a single 
Cu02 plane and fj,, v refer to the direction of the bond 
i.e. x or y. D lil ,(2q) = Z)^(R) exp (— i2q.R) is the 
static pair propagator of Cooper pair fluctuations with 
-D^(Ri M - Rjv) = {tyin^jv) ■ The quantity / M (k, q) = 
cos[(fc M — qfi)a] is a form factor arising from the coupling 
between a tight-binding lattice electron and a nearest- 
neighbour bond pair. Because of the d-wave LRO de- 
scribed as 'Ne'el' order [12] of the 'planar spin' ip m in 
the bipartite bond-centre lattice, the standard sublatticc 
transformation (i.e. ip m — ¥ tp m = A m exp(i# m ) where 
0m = @m f° r a^-bonds and 9 rn = 6 m + it for y-bonds) 
implies D XX (R) = D yy (R) = -D xy (R) = -D yx (R) = 
D (R) where D(R) can be written as 

D(R m - R„) = {i> m ) (<) + D(R m - R„) . (6) 

Here D(R) is the fluctuation term. The LRO part 
(ip m ) = Arf leads to a d-wave Gor'kov like gap with 
A k = (A^/2)(cos k x a — cos k y a)\ the corresponding elec- 
tron self energy is E(k, iuj n ) — A^/(iuj n + £k)- In widely 
used phcnomcnological analyses [23] of ARPES data, this 
form is used above T c with lifetime effects, both diagonal 
and off-diagonal in particle number space, added to E, 
i.e. 

E(k,w) - -iT 1+ ^f\. r (7) 

Here Ti is single-particle scattering rate and r is as- 
sumed to originate due to finite life-time of preformed 
d-wave pairs [23]. 
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We propose here that as described above, the electrons 
move (above T c ) not in a pair field with ci-wave LRO 
which decays in time at a rate put in by hand, but in a 
nearly static pair field with growing correlation length £. 
We assume, as appears quite natural for a system with 
characteristic length scale £ that D(R) ~ exp(— R/£) 
for large R, while £ diverges at T c . This natural form for 
D(R) is found, for example, in the Bcrezinskii-Kostcrlitz- 
Thouless (BKT) theory [36, 38-40] for two dimensions 
and in a Ginzburg-Landau (GL) theory for all dimen- 
sions. As described in Appendix E, it so happens that 
for large correlation lengths, E(k, iui n ) is nearly the same 
as that for preformed d-wave symmetry pairs. Below T c , 
D(R) decays as a power law (i.e. D{R) ~ R~ v with 
t] > 0) due to order parameter phase or 'spin wave'-like 
fluctuations (see the discussion in Section III). We find 
here the consequences of these for the spectral function, 
i.e. 

A(k, ui) = - -Im [G(k, iui n -> w + iS)] , (8) 

7T 

measured in ARPES [4]. G(k, iuj n ) is obtained from 
Eq.(4) with the self energy calculated from Eq.(5) using 
the aforementioned forms of the pair propagator D(R) 
[or -D(R)] for temperatures above and below T c . 

III. ELECTRON SELF ENERGY AND 
PHENOMENOLOGICAL INPUTS 

The tight binding lattice Hamiltonian T-L is given by 
Eq.(2). Near T c the time dependence of ipij can be ne- 
glected (except when T c itself is close to zero where time 
dependence of ipij has significant effects) so that this be- 
comes 

W = - E Ujalau, - J2 ^ b h + M ■ (9) 

ij,o <ij> 

Eq.(9) describes electrons moving in a static but (in gen- 
eral) spatially fluctuating pair potential ipij whose cor- 
relation length (in our case) diverges as T —> T c . The 
effect of the long distance fluctuations on the electrons 
is captured by the self energy shown in Fig. lb and alge- 
braically described by Eq.(5), with the significant fluctu- 
ation wave vector q ~ £ _1 << a -1 . In a regime where 
the fluctuations in the real pair magnitude A m are short 
ranged, 

D{R) ~ (A(R))(A(0))(e i[e " (R) - e " (o)] ) = A 2 D(R) 

(10) 

for large R (R >> a) (as evident from Eq.(6), D(R) ~ 
A 2 {D(R) - | < e iS ^ > | 2 )). This decoupling between 
magnitude and long distance phase correlations is accu- 
rate for x £ x pt, a manifestation of which is the sep- 
aration between T* and T c . Our calculations based on 
Eq.(10) are therefore reliable in this doping range. We 



use the general form [36] D(R) = (Ai?)" T ' cxp (-R/£) 
(with A ~ a^ 1 ) in Eq.(5), expanding £k-2q and / M (k, q) 
in powers of q for qa << 1. The self energy S(k, iuj n ) 
(Eq.(5)) is (see Appendix B) then 

f Av k J (vk/£ - i sgn(uj n )(iuj n + ^ k )) 

(11) 

In this equation, T is the well known gamma function 
and Vk = ^7jj|r (with = |vk|) is the velocity (ex- 
pressed in units of energy) obtainable from the energy 
dispersion £k- There is a relatively small particle-hole 
asymmetric part that adds to A 2 , (in the numerator of 
Eq.(ll)) a term proportional to (uo n + £k), and is ig- 
nored henceforth. The above self energy does not affect 
the nodal quasiparticles owing to the k dependence of 
Ak = (A/2) (cos k x a — cos k y a) (of which the second term 
arises from the form factor appearing in Eq.(5)). Here 
it is very important to mention that the quantity Ak 
should not be confused with the spectral gap Ak that 
we are going to define below from the position of peak 
(as a function of w at a fixed k) of the spectral function 
A(k,u). The spectral function A(k, to) (Eq.(8)) is calcu- 
lated using the self energy of Eq.(ll), and the consequent 
gap (Ak) in the electronic spectra is in general different 
(except at T = 0), due to pair fluctuations, from the in- 
put Ak (or A that appears in Eq.(10)). An additional 
term describing the coupling of electrons to short range 
phase fluctuations (present in our theory, but not dis- 
cussed here) is probably relevant for nodal quasiparticles. 
The self energy of Eq.(ll) has a nonvanishing imaginary 
part as u> — > because of the decay of an electron into 
one of opposite momentum and small q Cooper pair fluc- 
tuations (see Fig.l b). Thus the electronic system is a 
non Fermi liquid. 

In order to plot the spectral peaks, e.g. in Fig. 2a, we 
use a small lifetime broadening (~ 3 meV, which is much 
less than the typical instrumental resolution of ~ 10 meV 
of ARPES [4, 5]). The nodal peaks, which would be 6- 
functions in its absence, can then be detected clearly. 

Above T c , we use in Eq.(ll) the BKT form 

f(aj,T) ~ aexp[b'(x)/y/T/T c {x)-l] (12) 

and r\ = r/ c = T c /(2irp s (T c )) = 0.25, corresponding to 
the universal Nelson-Kosterlitz jump value [36], the BKT 
transition temperature (say Tbkt) has been taken to be 
actual observed [41] T c for the cuprates which have small 
interlayer coupling. The observed T c can be substan- 
tially higher (~ 5 — 10 K) than the underlying BKT 
transition temperature. This temperature Tbkt, though 
a fictitious transition temperature that can only be ob- 
served by switching off interlayer coupling between Cu02 
planes, is expected to control the correlation length above 
T c , away from a narrow 3D-XY critical region [42]. The 
use of experimentally measured T c in Eq.(12) as well as 
the precise value of the exponent b'(x) are not crucial for 
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the present calculations, although certain details, e.g. the 
T- dependence of the Fermi arc length close to T c (within 
~ 10 K), can be sensitive to these. For the 2D-XY model 
[43] b' ~ 1.6. No systematic estimate of either £(x,T) 
or b'(x) exists for the cuprate superconductors. Hence 
we work with b'{x) obtained using the phenomenologi- 
cal functional of Ref.[12] (see Appendix D for details). 
The Kostcrlitz RG equations [39] predict the asymptotic 
BKT form for £(T) to be valid only over a narrow criti- 
cal regime above T c . However, in practice, the BKT form 
is found [43] to fit the Monte Carlo data for £(T) for a 
2D-XY model very well over a rather large region of tem- 
perature above T c , till (£/&) — 2. In the same spirit we 
use the formula of Eq.(12) over a wide range. 

Below T c , we have used the self energy of Eq.(ll) which 
has been evaluated using the long-distance power-law 
form D(R) ~ R~ n appropriate for a quasi-LRO state in 
2D with r\ = (T/2irp s ), as well as using a D(R) in which a 
small intcrlayer coupling is incorporated through a sim- 
ple anisotropic 3D 'spin wave' approximation (see Ap- 
pendix C). Both approximations produce similar results 
for A(k,Lj). The quantity T](x,T) = T/(2irp s (x, T)), 
with p s (x,T) being the supcrfluid density, is an input 
to both of these and we use p s {x,T) estimated from 
penetration depth data[44] for Bi2212. In addition, the 
anisotropic 3D form for D(R) contains explicitly the in- 
tcrlayer coupling characterized by the c-axis supcrfluid 
density p% or the anisotropy ratio e = {p%j Ps) 1 ^ 2 (~ 10~ 2 
for Bi2212[45]). Close to but below T c , the experimental 
values for r\ exceed (see Fig. 5 a) by a large amount the 
upper limit of n for a pure 2D system i.e. ry c = 0.25. This 
fact demonstrates the essential role played by the inter- 
layer coupling in these systems (the observed LRO below 
T c is because of it!). Also, the 2D-BKT results that we 
use above T c , are not reliable, due to the importance of 
3D fluctuations for a narrow temperature regime close to 
T c such that £(T) > e~ 1 c (c is the intcrlayer spacing); 
we compare our results with ARPES data outside this 
range, i.e. outside the critical regime. 

Here, it should be mentioned that wc have neglected 
Coulomb interaction between Cooper pairs, which are 
charged objects, and associated quantum phase fluctua- 
tion effects, which can be specially prominent near hole 
densities (extreme overdoped and underdoped) where T c 
vanishes. The power law form used below T c in our calcu- 
lation arises due to longitudinal phase fluctuations and 
Coulomb interaction can push these up to plasma fre- 
quency, which is quite large (~ 1 eV for cuprates [46, 47]). 
In this case, the longitudinal phase fluctuations behave 
classically and are relevant only above a quantum to clas- 
sical crossover scale T c \. In the absence of dissipation T c \ 
has been estimated to be around 100 K, although dissi- 
pation can reduce T c \ to a much lower value, ~ 20 K, as 
has been inferred in Ref. [47]. In our calculation, below 
T c , the effects (e.g. the deviation of the k-space gap struc- 
ture from d-w&ve form) of classical phase fluctuations in 
electron spectral function are considerable only for much 
higher temperatures, typically close to T c where T > T c \, 



Hence we expect the calculated effects to be genuine and 
natural consequences for phase fluctuating superconduc- 
tors such as cuprates, and not to be changed much by 
quantum fluctuation effects except near the end points 
of T c (x) curve. 

The other phenomenological input, namely A(x,T), 
has been estimated from ARPES data[18, 48, 49] in the 
following way. (More details are given in Appendix E). 
Above T c , the self energy of Eq.(ll) can be written a 
simplified form for rj = (in the actual calculation, above 
T c we use r\ = r\ c = 0.25, the BKT jump value), i.e. 



Eflc,wl 



w + £k + i(«k/£) 



(13) 



The resulting peak positions of the spectral function 
A(k,uj) on the Fermi surface (i.e. £k = 0) are at A k 
which can be easily obtained as 



A k = 



At = A 



-(- 



\V2~t; 



for |A k | < 



otherwise 



V2£ 



(14a) 



(14b) 



The above expressions imply that, above T c , there is a 
gapless (i.e. spectral peak at zero energy) portion corre- 
sponding to the k values such that 



|A k (T)| 



< 



(15) 



As evident from this relation, the extent of this gapless 
portion along the Fermi surface is temperature dependent 
since both Ak and £ depend on temperature. This is 
the phenomena of 'Fermi arc' [17]; we discuss it in more 
detail in the next section (Section IV). This criterion 
(Eq.(15)) enables us to deduce a temperature, T an , at 

which the antinodal gap, A an = \J Al n - (v an /(\/2£)) 2 , 

gets completely filled in, i.e. A an (T an ) = so that the tip 
of the Fermi arc reaches the antinodal Fermi momentum 
k an (here A an = A kan and u an = Vk an )- This implies a 
self-consistent condition for T an 



A(x,T an ) ~ | A 

an 



v an (x) 



V2£(x,T an ) 



(16) 



With the aid of the above, we estimate the phenomeno- 
logical input A(x,T) at the ARPES-measured pseudo- 
gap temperature T* [48, 49] by identifying it with T an , 
i.e. A(T*) ~ v an /(V2£(T*)) (the right hand side of this 
relation is already known as we have estimated £ from 
Eq.(12) and is obtained from the band-structure £k)- 
Also at T = 0, A(0) can be deduced from the zero tem- 
perature antinodal gap A an (0) measured in ARPES [49] , 
as there A an (0) = (A(0)/2)(cos(fc alli? a)-cos(fc an; j / )). We 
use a simple interpolation formula, A(x, T) = A(x, 0)[1 — 
sinh(a(a;)r)] and utilize the knowledge of A(T) at these 
two temperatures, namely at T = and at T = T* , 
to determine the parameter oi(x) (see Appendix E). 
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This chosen form for A(x,T) implies that A an (T), for 
a particular x < x opt , varies substantially with T only 
near T = T* ~ T an . This choice qualitatively mim- 
ics the experimentally observed temperature-dependence 
[17, 22, 48] of the antinodal gap (see Appendix E). 
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FIG. 2: The Fermi surface and the spectral density 
along it. a, Fermi surface (FS) in the first quadrant of the 
Brillouin zone for x ~ 0.11, a typical value. The angular 
position <j> on the Fermi surface is shown. The antinode cor- 
responds to <j> = 0° and node to (j> = 45° . b-f, Single-particle 
spectral density ^(k,^;) at 15 equally spaced points (red cir- 
cles) on the Fermi surface in a shown vertically shifted from 
node (bottom) to antinode (top), as a function of tempera- 
ture at x ~ 0.11. Blue and red lines correspond to Fermi arc 
and gapped portions of FS, respectively. 



IV. RESULTS 

The results for spectral density are described below 
and are compared with ARPES experiments for u> 100 
meV in I^S^CaC^Os+s (Bi2212) (with maximum T c , 
ymax ^ 92 K occurring at a; op t ~ 0.16) across the un- 
derlying Fermi surface (FS) (see Appendix A). As men- 
tioned in the preceding section, pairing fluctuation re- 
lated quantities such as r](x 7 T) [77 = (T/(2np s )) below 
T c where p s (x,T) is the measured superfluid density], 
£(x,T) as well as the average local gap [12] A(x,T) = 
(A m ) act as inputs to our calculation. The GL-like the- 
ory of Rcf.[12] provides a unified approach in which all 
the above quantities and many other results emerge from 
a single assumed free energy functional; this also leads 
to results [12] similar to those described below. As re- 
marked earlier, our results are not reliable for extreme 
undcrdoped and overdoped regions where T c and 
consequently quantum phase fluctuation effects can be 
crucial [12, 35]. 

Our results for the (particle-hole symmetrized) single 



electron spectral density A(k, uj) are shown in Fig.2b-f. 
From the spectral density curves (Fig.2b-f) for an un- 
dcrdoped system (x ~ 0.11 and T c ~ 67 K) it is clear 
that above T c , for some region of the Fermi surface cen- 
tered around the nodal point (0 = or k x = k y ), A(k., lu) 
peaks at lu = 0. The extent in k space over the Fermi 
surface of this (ungapped quasiparticle excitation) part, 
normally identified with the 'Fermi arc' [17], depends on 
temperature. For any k outside this region, the spec- 
tral function has a peak at a nonzero energy; the value 
of the peak position is identified conventionally with the 
energy gap Ak- The Fermi arc emerges above T c due to 
the existence of a finite correlation length in the system, 
as explained below. 

As discussed in the preceding section (Section III), the 
self energy expression (Eq.(ll)) obtained by us implies 
that the gapped portion terminates (and the Fermi arc 
begins) when |Ak| < fk/£- One can understand this re- 
lation from the fact that a quasiparticle with momentum 
k near the Fermi energy moving with a velocity Vk in 
a spatially fluctuating superconductor having correlated 
patches of length £ undergoes scattering over a time in- 
terval t £ £/vk- The energy of such a quasiparticle is 
~ Ak so that from the uncertainty relation rAk ^ 1 
or (£/ w k)Ak £ 1. Above T c , even in the region of the 
Fermi surface where Ak 7^ 0, there is (both in our calcu- 
lations and in experiments) substantial spectral density 
at zero excitation energy. The simultaneous presence of a 
nonzero Ak and finite spectral density at zero frequency 
is the operational definition of pscudogap. Below T c , the 
calculated spectral density has two symmetrical peaks at 
nonzero Ak for all k except at the node, where the peak 
is at zero energy. 

Fig. 3a gives an overall pictorial view of the develop- 
ment of the Fermi arc as a function of x and T above 
T c . In Fig. 3b, we plot the arc length versus the reduced 
temperature, namely (T/T an ), where T an is the temper- 
ature at which antinodal pscudogap fills up. The nearly 
straight line obtained by us over a large range of doping 
(0.07 < x < 0.15) is seen to compare well with ARPES 
data over the same range [18]. Qualitatively, since the k 
points on the Fermi arc satisfy the condition | Ak| ^ Vk/£, 
this condition is met for a larger k as £(T) decreases on 
increasing T. Thus the arc length increases with increas- 
ing T, and at T an {x) ~ T*(x) it is as large as the entire 
Fermi surface. Conversely, it shrinks to zero at T c and 
below because £ — > 00 there. 

Fig. 3c shows T an (x) together with T c (x) [41] and ex- 
perimental data for T*(x) [49] for Bi2212. As men- 
tioned above, T an is taken to coincide with the tem- 
perature at which the Fermi arc takes up the full FS, 
i.e. |A an (z,T an )| = v an /(V2Z(x,T an )). Hence T an ( x) is 
decided by A(x,T) ~ |A an (x,T)| and £(x,T), both of 
which are phenomenological inputs to our calculations 
(v an and k an are obtained from the energy dispersion 
£k). For £(x,T) we use the BKT form of Eq.(12). We 
have already discussed the choice of the exponent b'(x) 
appearing in Eq.(12) in Section III and we take a to be 
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FIG. 3: Spectral properties above T c . a, Colourmap of A{k,ui = 0) over the first quadrant of the Brillouin zone for 
two values of x at different temperatures. The Fermi arc is easily picked out visually. The colourbar indicates the value of 
A(k,uj = 0) (in units of 3 x 10" 3 meV" 1 ). b, The arc length vs. T/T an curve excluding the temperature region within ~ 15 
K of T c (see Appendix D), averaged over the entire doping range (x < x opt ). On the y-axis, 0% is the node and 100% is the 
antinode. The experimental data are from Ref. [18]. The vertical error bars in the theoretical points indicate the variation of 
the arc length at different x for a given T/T a _ n (x). c, The antinodal pseudogap filling temperature T an as a function of x is 
compared with the data for the pseudogap temperature T* from ARPES [48, 49]. The parabolic T c (x) curve [41] used in our 
calculation is also shown. 



the square planar Cu lattice constant, thus completely 
fixing £(x,T) for a given x and T. The main features of 
our results are not crucially dependent either on the spe- 
cific form of £(2;, T) used here or on the particular chosen 
value of b'(x) (see Appendix D). The other input A(a;, T) 
(Appendix E) is fixed by enforcing T an , obtained via the 
relation A(x,T an ) ~ v an /(\/2£(x, T an )), to be close to 
the ARPES data [49] for the antinodal pseudogap filling 
temperature T*(x) i.e. the pseudogap temperature scale 
measured in ARPES. Thus T an (x) is preordained to be 
nearly same as T*{x) in our calculation, Fig. 3c manifests 
this fact. 

Some other spectral features observed as T decreases 
through T c are shown in Fig. 4. Each of the left figure 
panels (Fig. 4a, c) exhibits a particular spectral prop- 
erty (see below) over the entire FS (see Fig. 2a), for a 
typical x ~ 0.11 at several temperatures. In the respec- 
tive right panels (Fig. 4b, d) we show the same spec- 
tral quantities specifically at the antinode (4> = 0) as 
a function of temperature for various x. The pseu- 
dogap is widely characterized [18, 48] by the function 
L(4>) = 1 - A((f),0)/A(<p, A) that quantifies the loss of 
spectral weight A(cj), lu) for energy lu = ±A(0) at angu- 
lar position (f> (Fig. 2a) on the FS. If there is a nonzero 
gap in the excitation spectrum, and there are no zero en- 
ergy excitations (as in a BCS superconductor at T = 0) 
L((p) = 1; if there is no gap (real or pseudo) L(<fi) = 0. 
Intermediate values (i.e. < L(4>) < 1) imply a pseudo- 
gap. The evolution of L((p) as T decreases through T c 
describes how the pseudogap develops into a real gap. 
This is shown in Fig. 4a. It agrees qualitatively with 
that obtained from ARPES [48]. Experiments show a 



nonzero A(k, lu = 0) (or L(<fr) < 1) below T c as well, un- 
like in our theory where below T c , L((f>) = 1, except at 
(f> = 45° where it is zero. The difference is presumably 
due to neglect of some causes of intrinsic spectral line- 
shape broadening, not included in our model, and finite 
instrumental resolution of ARPES. Above T c as well, the 
calculated L((f>) in the pseudogapped portion is less (by 
about 20%) than the experimental values, probably due 
to such effects. In Fig. 4a, we have (arbitrarily) multi- 
plied our results by a related factor (~ 1.2) for compari- 
son with experiment. The spectral loss at the antinodal 
point i.e. L(<p = 0) = L an is exhibited in Fig. 4b for sev- 
eral values of function of T. L an decreases almost 
linearly with temperature in line with observations [48] 
and vanishes at T = T an . 

Another aspect of the spectral function A(k, u) is 
shown in Fig. 4c, in which the full width at half maximum 
(FWHM) of the peak is plotted as a function of <p, again 
for x ~ 0.11. We notice a maximum in it which coincides 
with the 4> location of the Fermi arc tip. Similar features 
can be deduced from the ARPES spectra[48]. We also 
find in Fig.4d that while the FWHM of the antinodal 
peak increases substantially with T below T c in the un- 
derdoped case, it does not change much for the overdoped 
cuprate, slightly beyond optimal doping. This is what is 
observed[19]. Again, the experimental widths are seen 
to be around 20 meV larger than the calculated ones, 
perhaps due to the same neglect of some intrinsic and 
extrinsic quasiparticle lifetime effects mentioned earlier. 

Possibly the most widely explored property of a super- 
conductor is the energy gap Ak or A(0). The apparent 
deviation ('bending') of the gap below T c , inferred from 
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FIG. 4: Spectral properties both above and below T c . 

a, Loss of low energy spectral weight L(<f>) (see text for defi- 
nition) for three temperatures at x ~ 0.11. The symbols cor- 
respond to the experimental data [48] at same temperatures 
for an underdoped (UD) sample with T c = 67 K, T* ~ 150 
K. b, Filling of the antinodal pseudogap (quantified by L a n) 
with temperature for three values of x. c, Width (FWHM) 
of the spectral peaks along the FS for x ~ 0.11. d, Variation 
of the width of the antinodal peak with temperature (< T an ) 
for an underdoped (x ~ 0.12) and an overdoped (x ~ 0.19) 
sample, compared with antinodal widths at the same x values 
extracted from ARPES data in Refs.[19, 20]. Vertical dashed 
lines indicate the corresponding T c 's. 
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FIG. 5: Spectral properties below T c . a, Deviation of 
Ak from the d-wave form, \ cos k x a — cos k y a\, below T c for 
x = 0.16 is compared with Ak obtained in ARPES [22]. At 
low temperature (T = 10 K) , Ak follows the canonical d-wave 
form (a straight line), but bends away considerably from it 
close to T c . For T < T c , we have shown Ak obtained for both 
the true LRO phase (i.e. D(R — s> oo) = \{ipm)\ 2 0) and 
the quasi-LRO or QLRO phase (i.e. D(R) ~ IT" 0) 
(see Methods). The quantity r\ oc T/p s or the superfluid den- 
sity p a {x,T) has been estimated from the penetration depth 
data of Ref. [44]. b, The antinodal peak height A? n vs. T 
tracks p 3 (T)/T vs. T. Inset At low temperatures (T — 6 
K) A?, n (x) closely follows p s {x), especially in the underdoped 
side {x < 0.16). 



ARPES [21, 22], from the d-wave form has been the sub- 
ject of a great deal of current interest, leading to specula- 
tion that there are two gaps in high-T c superconductors 
[50]. We have obtained the gap Ak from the peak posi- 
tion of the calculated spectral function A(k, uj) for a fixed 
k as mentioned earlier and conclude from the results (see 
Fig. 5a) that the 'bending' is due to the coupling of the 
electron to thermal phase fluctuations ('spin waves') be- 
low T c . As expected from such an origin, it is large close 
to T c , and small as T — > 0. These results confirm that 
there is only one gap, but that to 'uncover' it, effects of 
coupling to pair fluctuations ('spin waves') have to be 
included. 

Below T c , because of the LRO in the Cooper pair am- 
plitude, the electrons move in a lattice periodic pair po- 
tential while occasionally getting scattered from the ther- 
mal 'spin wave' fluctuations around the ordered state, so 
that the eigenstatcs have coherent quasiparticlc features. 
One consequence, shown in Fig. 5b, is that the height 
(A^ n ) of the coherent antinodal peak at uj = ±A an follows 
closely the T-dependence of r\~ x oc (p s /T). A similar em- 
pirical correlation has been reported in Refs.[19, 20]. At a 
given temperature, A? n is proportional to p s {x 1 T). This 
is shown in the inset of Fig. 5b for a particular value of 
T as a function of x. 



V. DISCUSSION 

Our work is a comprehensive exploration of the con- 
sequence of the (inevitable) coupling between electrons 
and Cooper pair fluctuations constituted of the same 
electrons in a two dimensional model for the cupratcs 
with a square lattice. Underlying this is a picture of the 
superconducting transition as a continuous approach to 
long range order characterized by nonzero phase stiffness 
so that the pair (phase) correlation length £ diverges at 
T c . (One way this can happen, which is the mechanism 
discussed in an earlier paper [12] by us in a Ginzburg- 
Landau-like theory of superconductivity in the cuprates, 
is the following. Nearest neighbor spin singlet pairs are 
the basic degree of freedom, and are preformed well above 
T c for optimal and suboptimal hole density. Their near- 
est neighbor phase dependent interaction of 'AF' type 
leads to d-wavc symmetry LRO below T c ). 

We have obtained here explicitly the self energy 
E(k, iu> n ) of an electron moving in the field of Cooper 
pairs whose d-wave symmetry correlation length £ di- 
verges as T — » T+. This is done very generally, e.g. as- 
suming that the large distance (R >> a) phase corre- 
lator function is of the form R~ v exp (— where ij 
is the anomalous dimension. In actual calculations we 
use the specific forms (and values where appropriate) for 
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T] and £ from the BKT theory [38-40] without making 
any further assumptions regarding rj and £. Since the 
BKT theory predicts somewhat lower T c compared to 
experiment (T c — Tbkt ~ 5 — 10 K) and since the effect 
of small interlayer coupling which leads to the observed 
three dimensional (XY) critical behavior very close to T c 
is neglected, we compare our results for spectral func- 
tion A(k, u) with ARPES experiments outside a regime 
5T ~ 10 — 15 K from T c . In this paper, we concentrate 
on ARPES results because ARPES represents real space 
averaged measurements. We do not address STM (Scan- 
ning Tunnelling Microscopy) results because the latter is 
a spatially local measurement sensitive to local inhomo- 
geneities. 

The lowest order vertex correction to the self energy 
shown in Fig. la vanishes identically above T c . This is 
because the vertex correction involving the Cooper pair 
correlation function (ip m ^) requires the anomalous elec- 
tron propagator to be non zero. The latter is identically 
zero since there is no long range superconducting or- 
der. The part of the vertex correction involving (ipm^n) 
(or (ipm^n)) is zero because such functions vanish due 
to 'isotropy' or 'XY-symmctry' in the space of two- 
component 'spin' ip m = (A m cos# m , A m shi9 m ). Below 
T c , the correction can be calculated exactly as in the clas- 
sic work of Migdal [51]. The internal phonon propagator 
of the Migdal calculation is replaced here by the collective 
phase fluctuation or spin wave or Goldstone mode prop- 
agator, so that the correction is ~ (T c /ep) ~ 10~ 2 << 1. 
However, higher order vertex corrections are nonzero. 
These have not been calculated. The internal 'bosonic' 
or pair fluctuation propagator is the true one, not the 
bare one, since it is described in terms of the observed 
transition temperature T c which fully includes dressing 
or renormalization effects . The approximation made is 
that the bare fermion (electron/hole) propagator is used 
in calculating the self energy instead of the dressed one. 
Experience with the Eliashberg approximation in phonon 
induced Cooper pair theory (BCS) as well as calculations 
for coupled electron phonon systems in the normal state 
show that this is reasonable. We would like to point out 
that the role of the 'boson' here is not to form Cooper 
pairs (as in the celebrated BCS theory); the spin singlet 
nearest neighbour Cooper pairs are already formed while 
the interaction between them leads to the emergence of 
d wave symmetry superconductivity (long range order 
or LRO) below T c as a collective effect. The 'boson' of 
Fig. la is the correlation at large distance between d-wave 
symmetry or collective Cooper pair fluctuations, whose 
length scale diverges at T c . 

We do not attempt to explain the broad part of the 
one electron Green function [4. 30, 49], specially promi- 
nent away from the node, both above and below T c (this 
'incoherent' part may well be a strong correlation ef- 
fect). We assume however that at low excitation ener- 
gies there is a 'quasiparticlc' part in the 'bare' one elec- 
tron Green function, with weight Zk (the corresponding 
Green function, that we calculate here, only gets scaled 



by Zk, i.e. G(k, iw n ) ZkG(k, iuj n ), so that the area un- 
der the quasiparticlc part of the spectral function A(k, u) 
is Zk [30]). In strongly correlated system Zk is usually k- 
independent [30]. Recent QMC simulations [52, 53] of 
the one band Hubbard model in two dimensions (widely 
believed to be appropriate for the cuprates) shows that for 
strong correlations, there is indeed a quasiparticle part 
to the one electron Greens function, and that the quasi- 
particle residue Zk is weakly k dependent. Presumably, 
such a z can be absorbed into a definition of renormal- 
ized A 2 , i.e. (A 2 ) = (A 2 z) while evaluating G(k, u>) 

using the Dyson equation (Eq.(4)) and the self energy ex- 
pression of Eq.(5), as G°(k,u) -> zG°(k, uj) there. The 
electronic self energy we calculate as a result of coupling 
to pair fluctuations depends strongly on k; in particular 
it vanishes at the nodal point. 

There have been a large number of calculations of the 
electron spectral density in cuprates under more specific 
assumptions and addressing particular experimental find- 
ings, e.g. the Fermi arc or the pseudogap. We mention 
below some that we are aware of. None, as far as we 
know, compare theoretical results with those of actual 
experiments explicitly, over this wide a range of temper- 
ature, doping, experiments and phenomena or uses this 
picture of the superconducting phase transition. 

Many earlier calculations of the effect of phase fluctu- 
ations (e.g. Franz and Millis [24], Berg and Altman [25]) 
assume that these are connected specifically with super- 
currents surrounding thermal vortices. While we use the 
BKT model, which describes the superconducting tran- 
sition as due to vortex-antivortex unbinding, and the as- 
sociated correlation length £ (Eq.(12)), our form for the 
self energy (Eq.(ll)) is independent of the actual mech- 
anism for the phase fluctuation spectrum or the way £ 
diverges approaching T c , depending only on the gener- 
ally valid form for the spatial dependence of of the cor- 
relation function of the phase fluctuations. Some other 
recent calculations (Micklitz and Norman [55], Senthil 
and Lee [54]) either use different forms for D(R) or as- 
sumes implicitly that n = and £ = oo. The latter 
implies preformed d-wave symmetry pairs which is al- 
most the prevailing belief in the field. In our approach, 
eZ-wavc symmetry long range order develops at T c ; above 
it this begins to emerge as a collective phenomenon as 
£ — >• oo. A recent calculation by Tsvelik and Essler [56] 
is similar in spirit to ours, but in our opinion it is too 
strongly tied to the anisotropic Kondo limit. Perhaps 
the closest is the recent numerical calculations by Li and 
coworkers, who have analyzed the electron self energy £ 
via Monte Carlo sampling of 2D-XY model [57] as well 
as the Hubbard-Stratonovich transformed version [58] of 
the pair attraction term of Eq.(2) using an involved nu- 
merical technique. Their results are similar overall to 
some of those described here, e.g. in Fig. 3a. 

Our approach for calculating the low energy excita- 
tion spectrum of electrons coupled to quasi static long- 
wavelength fluctuations implies that an electron decays 
into another with nearly opposite momentum and a 
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Cooper pair fluctuations with small momentum 2q. This 
is taken to be the only contribution to quasiparticle decay 
near the Fermi energy. We do not obtain or emphasize 
finite Cooper pair life-time effects, generally called 'dy- 
namic' effects in the literature, and mostly included em- 
pirically to fit ARPES data. We also do not use results 
from BCS theory which being a mean field theory, can 
not, wc believe, describe accurately the region above T c 
(and say, below T*) since here local Cooper pairs (near- 
est neighbor spin singlets) exist, but condense collectively 
into a state of d-wave symmetry LRO only at T c . 

As mentioned earlier, a variety of exotic mechanisms, 
such as stripes [16], d-density wave [13] and time rever- 
sal symmetry breaking orbital current state [14], have 
been suggested and invoked to describe the phenomena 
observed in the cuprates (especially above T c ). Our work 
provides an explanation in terms of a simple, inevitable 
process. We believe that the results, and their agreement 
with experiments, strongly support our mechanism [12] 
of nearest neighbor Cooper pairs, and long-range d-wavc 
symmetry order emerging as a collective effect arising 
from short-range interaction between these pairs. This 
probably points to the way in which, we believe, high-T c 
superconductivity will be understood. 

We do not include here the effect of coupling to ad- 
ditional short-range fluctuations which we find to be 
present. These may have a decisive effect on nodal 
quasiparticle life-time as well as on linear resistivity in 
the strange metal phase. We have also not emphasized 
particle-hole symmetry breaking terms in A(k, ui) which 
we have calculated. Our approach implies that fluctua- 
tion effects, very small in conventional superconductors, 
are sizable in the cuprates and are present over a wide 
temperature range. 
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Appendix A: Determination of the Fermi surface 

The chemical potential /i is computed by setting 
no^Ek/H^w) = (!-*) (herc/H = l/(e^ + 
1) is the Fermi function) and the underlying Fermi sur- 
face is defined from the locus of £k = 0. We have im- 
plemented other FS criteria in our calculation e.g. the 
locus of the minimum of At [18] as well as the locus 
of the maximum of A(k, tu = 0). The main features of 
the results (Figs. 2-5) remain unaltered for different FS 
criteria. 



Appendix B: Evaluation of self energy 

When the pair correlator D(R) in Eq.(6) is sufficiently 
long-ranged (below T c at all temperatures, and above 
T c for temperatures such that £(T) >> a), its Fourier 
transform D(2q) = £>(R) cxp (— i2q.R) is sharply 
peaked around q = and one can expand £k-2q and 
//j(k, q) in Eq.(5) in powers of q, so that £(k, iu> n ) can 
be approximated as, 

£(k, iu> n ) ~ 7?.(cos k x a — cos k y a) 2 + . . . , (Bl) 

where K = A" 1 £ q D(2q)/(iw„ + £k - 2v k .q). In 
Eq.(Bl), we have not shown explicitly particle-hole non- 
symmetric contributions to S(k, iuj n ) for which we have 
obtained closed form expressions. This can be recast us- 
ing standard mathematical identities into forms suitable 
for analytical and numerical calculations, i.e. 

U = - l sm ^ n) /„°° dse is ^M^+^>D{sv k ), (B2) 

The above expression utilizes only the real-space pair cor- 
relator D(R) (Eq.(6)). This procedure is very different 
from the way self energy is generally calculated using mo- 
mentum space form of D(R) (see e.g. Refs.[54, 55]). For 
instance, using D(R) ~ A 2 (T)D(R) (Eq.(10)) along with 
D(R) = (AR)-^cxp(-R/£_) (A = e - Q - lle 2V2^/a is re- 
lated to upper wave- vector cutoff [36]) results in the self 
energy expression of Eq.(ll). As evident, the integration 
variable s in Eq.(B2) has dimension of time or inverse en- 
ergy and hence the presence of D(sv-k) ~ exp (— si>k/£) 
in Eq.(ll) cuts off the contributions for large timescales 
corresponding to s » £/«k- Below T c , D(R), calcu- 
lated from a spin wave approximation incorporating fi- 
nite interlayer coupling (see below), is used in Eq.(B2) 
to obtain S(k, iu> n ) by numerical integration. This also 
demonstrates the usefulness of Eq.(B2) for evaluating self 
energy using arbitrary form of the real-space pair corre- 
lator as it amounts to carrying out a one dimensional 
integral, cither analytically or numerically. 



Appendix C: Pair correlator f (R) below T c in the 
presence of non-zero interlayer coupling 

Below T c , we estimate the phase correlator D(R) = 
^ e j(e(R)-0(o))^ Q £ Eq.(6) using the well known harmonic 
spin wave approximation i.e. the elastic free energy func- 
tional for an anisotropic 3D system, 

•Fsw = ^ J dr{(d x <p) 2 + (d y cp) 2 } 

+ J dr(d z <pf, (CI) 

where d x = and c is the interlayer spacing. As men- 
tioned earlier, p s = p® b and p c s are the a6-plane and c-axis 
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using the anisotropic 3D calculation is only accurate at 
very low temperature for large anisotropy or small e. In 
our calculation, this fact is manifestly observed e.g. in 
the variation of the antinodal gap A an with temperature 
for T < T c (see Fig. 8b). The planar supcrfluid density 
p s {T) is estimated from the measured afo-plane penetra- 
tion depth [44] for Bi2212 using p s = $%c/(l6n 3 \ 2 ab ); $ 
is the fundamental flux quantum and we take c = 15 A as 
suitable for Bi2212. For the anisotropy ratio e (assumed 
to be temperature independent), we use the empirical 
formula [45] e(x) = e(x op t)(a; - x u )/(x opt - x u ), where 
a^opt = 0.16 is the optimal hole doping and x u = 0.047 is 
the end point of the T c (x) dome in the underdoped side, 
with e(xopt) — 1/133. 



FIG. 6: Phase correlator D(R) below T c . D(R) (solid 
lines) estimated using Eq.(C3) for x — 0.124 at four different 
temperatures. Also shown are the corresponding 2D correla- 
tor D2o{R) (dashed lines). The vertical dashed line indicates 
the length scale below which D(R) ~ D2u(R)- 



supcrfluid density, respectively. The phase correlator can 
be written [36] as £>{R) = exp(—g(R)), with 



g(R) 



Tc 



7r / a 



dq x dq y r /c dq z 1 - e^ R 



Ps{T) J-Tt/a (27T) 2 J. %/c 27T ql + ql + eHcq z /a)^ 

(C2) 

Here e = (p c s / Ps) 1 ^ 2 and R = (R x ,R y ,0) is in the xy- 
plane. Finally one obtains, 

D(R) 
= exp 




Here, Jo is the ordinary Bessel function, and c/lro = 
2^[A7r-2Acot~ 1 (£7r/A)+e7rhi(l + (A/7re) 2 )] with A = 
and rj = T/(2irp s ). The phase correlator D(R) is 
shown in Fig. 6a for x = 0.124 at different temperatures. 
£>(R) starts from D(0) = 1 and decreases to a constant 
value -Dlrq = exp(— <7lro) as i? — ^ oo. The fluctu- 
ation part .D(R) (Eq.(6)) can be obtained as D(R) ~ 
A 2 (T)(D(R) — -Dlro)- In Fig. 6a, we also show the corre- 
sponding 2D phase correlator (Z?2d(-R) = (Ai?) - '') which 
is appropriate for a QLRO phase as £>2d(^ — > oo) — > 0; 
D-2T){R) is also obtained from the spin wave approxima- 
tion for a 2D system[36]. For a < R < e^a, D{R) ~ 
D2d{R) i-e. the system effectively behaves as two di- 
mensional. Both the spin wave approximations, the 
anisotropic 3D and pure 2D, are expected to be quanti- 
tatively correct only at low temperatures. Especially the 
estimation of the LRO part of D(R) i.e. -Dlro obtained 



Appendix D: Correlation length £(x,T) above T c 

As mentioned in Section III, £(x, T) is esti- 
mated using the BKT form [36] i.e. T) ~ 
aexp(b'(x)/y / T/T c (x) — 1). To obtain the actual cor- 
relation length, one needs the value of the dimension- 
less quantity b'{x), which is in the exponent. This is 
a nonuniversal number of order 1. For thin supercon- 
ducting films its value ranges [42] from 1 to 4. Unfor- 
tunately for cuprates no such estimate exists for b' (pre- 
sumably doping and material dependent). We have esti- 
mated b' using a GL-like functional in Ref.12. As shown 
there, with appropriate choices of a few parameters as 
input in the GL-like functional, it can be used to explain 
and reconcile a wide variety of experimental data for the 
cuprates e.g. those for T c (x), p s {x,T), and the contri- 
bution of pair degrees of freedom to the electronic spe- 
cific heat. The quantity b'(x) has been estimated there 
in the following manner. In two dimensions, the sys- 
tem described by the model undergoes a BKT transition 
at a temperature Tbkt(^)- So, the correlation length 
£(x, T) is expected to follow the aforementioned BKT 
form above T c (a;), where T c (x) is identified with Tbkt(x). 
In such a scenario, as predicted by BKT theory [36], the 
supcrfluid density below T c (x) is given by the formula 
[36, 59] p s (x,T) =p s (T c (x))[l + b(x)y/l-T/T c (x)] near 
the transition. The BKT theory further predicts a uni- 
versal relation[59] between b' and b, namely bb' = tt/2. 
Keeping this fact in mind, we calculate p s (x,T) below 
T c from a Monte Carlo simulation of the GL modcl[12] 
and fit the results to the above mentioned BKT form for 
p s (x,T) to estimate b(x) and thence b'(x). The results 
for b'(x) are shown in the inset of Fig. 7a). 

We show £(T), estimated using b'(x) obtained in this 
manner, in Fig. 7a for four different x values. These re- 
sults for T) have been used in Eq. (1 1) to calculate the 
spectral properties above T c {x), as reported in Figs. 2-5 
of the main paper. We have also used the constant value 
b'(x) = 1.59, expected for a 2D-XY model[43], to repeat 
the same calculations. The results in the two cases only 
differ in details and their main features over the entire 
x range are found to be robust with respect to different 
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FIG. 7: Correlation length £(x,T) and Fermi arc length, a, The correlation length £, which controls the approximate 
linear temperature dependence of arc length with temperature in Fig. 3 b, is shown as a function of dimensionless temperature 
T/Tan{x) for four x values. Inset: b'(x) estimated from the GL-like functional of Ref.12. The horizontal dashed line corresponds 
to b' — 1.59, as estimated for a 2D-XY model [43]. b, The arc length vs. T/T ali (x) curves (the average curve for the doping 
range 0.07 < x < 0.15) for two different choices of b' , namely 6' = 1.59 and an ^-dependent b' estimated using the functional 
of Ref.12. The rapid increase of arc length near T/T an = 1 is related with the rapid decrease of (A m ) = A(T) near T an (see 
Fig. 8). c, The arc length as a function of T/T an for x ~ 0.11. ARPES results[18] for the arc length for the same doping at a 
few values of T/T an are also shown. 



choices of b'(x) ~ 1 (sec Fig. 7b). These temperature de- 
pendence of £ _1 (along with (A m ) = A(T), see below) 
governs the behaviour of the arc length as a function of T 
(see Fig. 7c). Starting from zero at T c , the arc length rises 
steeply with T in a temperature range (~ 15 K) from T c 
(the actual value depending on x) to ultimately follow 
the mean curve shown in Fig. 7b. In our calculation this 
temperature range depends primarily on the details of 
the temperature dependence of £ close to T c and hence 
the range can be tuned by changing b' as well as the T c 
appearing in Eq.(12). 




100 200 300 50 100 150 

T (K) T (K) 



We use actual observed T c values for evaluating £ from 
Eq.(12). It is well known that the BKT theory result for 
T c (Tbkt) in principle can be substantially lower than 
the observed T c [42]. Also, we have neglected interlaycr 
coupling which affects £(T) and indeed determines the 
critical behavior very close to T c (It is 3D-XY [45] and 
this 3D critical regime is about a degree or two wide). 
Because of the two reasons, we compare our results for 
the Fermi arc length in Fig. 3 with experiment only be- 
yond a regime ST ~ 10 — 15 K from T c . Using our the- 
ory, £(x,T) itself can be estimated more accurately at a 
phenomenological level from a detailed comparison with 
experimental data for the arc length as a function of T, 
especially in a temperature region close to T c since A(T) 
is weakly temperature dependent in this regime. 



FIG. 8: Average local gap and antinodal gap. a, 

A(T) = (A m ) for four values of hole density. Vertical dashed 
lines indicate the respective antinodal pseudogap filling tem- 
peratures T an . b, A an (r) for the same hole densities. 



Appendix E: Fermi arc criterion, average local gap 
A(x,T) and antinodal gap A an (a-,T) 

As already discussed in Section III, the genesis of the 
Fermi arc above T c can be understood from the self 
energy of Eq.(ll). Neglecting the small rj in Eq.(ll), 
S(k, uS) = S(k, iu n — > uj + iS) can be written in a simple 
form depicted in Eq.(13). This is exactly same as the 
phenomenological form[23, 60] that has been used widely 
to analyze the ARPES data using a single particle scat- 
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tcring rate Ti and d-wave pairs with finite life-time Tq. 
In Eq.(13), To is replaced by v k /£ (the single particle 
scattering rate is zero for r) = 0), although our picture 
for its origin is completely different, as described in the 
main paper. 

For k lying on the Fermi surface (£ k = 0), the spec- 
tral density A(k, ui) has a peak at u = for k val- 
ues such that |A k (T)| < v k /(\/2£(T)) and two peaks 

at to = ±WA|(T) - (v k /(\/2£(T))) 2 otherwise. The 
former determines the extended region corresponding 
to the Fermi arc. The antinodal pseudogap A an ~ 

J A 2 n — (v an / (-\/2£)) 2 is completely filled in above a tem- 
perature T an , at which the end of the arc reaches the 
antinodal Fermi momentum [61] k an i.e. the criterion of 
Eq.(16) is satisfied. 

At this point, we would like to reemphasize that both 
A(x,T) and £(x, T) are purely phenomenological inputs 
to our theory. As mentioned in Appendix D, £(x,T) 
has been determined following fairly general considera- 
tions related to the superconducting transition at T c {x) 
(we use the empirical T c (x) curve [41] appropriate for 
Bi2212, see Fig. 3c). To choose A(x,T), we first demand 
T an (x) to be close to the pseudogap filling temperature 
T* obtained in ARPES experiments [48, 49]. Since v an (a;) 



(calculated from the tight-binding energy dispersion £ k , 
see Section III) and £(x, T an ) (Appendix D) are already 
fixed, Eq.(16) enables us to deduce the value of A(x, T a „). 
Also, at T ~ 0, A(x, 0) ~ A an (x, 0). Using this relation, 
we obtain the zero temperature antinodal gap A an (x,0), 
again from the ARPES data[49], to estimate the value 
of A(x,T) at T = 0. Once both A(x,0) and_A(x, T an ) 
are known, a simple interpolation formula, A(x, T) = 
A(x, 0)[1 — sinh (a(x)T)], is used to obtain the parame- 
ter a(x) so that the function A(x,T) is completely de- 
termined (Fig. 8). This interpolating form is chosen to 
qualitatively mimic the fact that A an (T) (~ A(T)) varies 
substantially with temperature only near T = T an , par- 
ticularly in the underdoped Bi2212, as reported in Ref.48. 
However, we note that very different temperature varia- 
tion of A an (T), i.e. A an (T) varying considerably with 
temperature below T an , has been reported in ARPES 
studies [62, 63] on other cuprates (e.g. Bi2201, La214). 
These observations suggest that the form of A(x, T) may 
differ substantially from one material to other. Since 
A(x,T) is a phenomenological input to our theory, such 
issues can not be resolved unless a satisfactory micro- 
scopic theory for cuprate superconductivity is developed 
such that A(x,T) is calculable starting from an appro- 
priate microscopic Hamiltonian. 
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